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ABSTRACT 

Aims. Analysis and interpretation of the Her X-1 X-ray light curve obtained with the ASM onbord RXTE over the period 1996 February to 
2004 September. 

Methods. We construct the averaged X-ray light curves by means of adding up light curves corresponding to different 35 day cycles. A 
numerical model is introduced to explain the properties of the averaged light curves. 

Results. We report that the features found previously in the averaged X-ray lightcurve are confirmed by the new RXTE/ASM data. In 
particular, anomalous dips and post-eclipse recoveries in two successive orbits in the short-on state are found to be prominent and stable details 
of the light curve. We argue that a change of the tilt of the accretion disk over the 35 d period is necessary to account for these features and 
show that our numerical model can explain such a behavior of the disk and reproduce the observed details of the light curve. 
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1. Introduction 

Her X-l/HZ Her is a close binary system with a 1.7 d or- 
bital pe riod consisting of an accretion powered 1.24 s X-ray 



pulsar dGiacconi et alj 1 1973c iTananbaum et alj 119721) and an 



optic al compan i on - a main sequence star of spectral type 
A/F dCramptonl 1 1 9741) . The X-ray flux of the source shows 
a ~35 d periodicity which is thought to be due to a counter- 
orbitally precessing tilted a ccretion disk around the neutron star 
dGerend & Bovntonlll976l ). 

The main feature of the X-ray light curve is the alternation 
between on (high X-ray flux) and off (low X-ray flux) states. 
The 35 d cycle contains two on states - the main-on (~7 orbital 
periods) and the short-on (~5 orbital cycles) - separated by ~4- 
5 orbital cycles. The maximum X-ray flux of the main-on is ~5 
times higher than that in the maximum of the short-on. The off 
states are usually explained by periodic obscuration of the X- 
ray source by the disk wh ich most likely has a twisted form 
dGerend & Bovntonlll976l) . 

According to the g enerally accepted model (see e.g. 



Schandl & Meyer 1994; Kusteretal. 2005, and references 



therein) the main-on state starts when the outer disk rim opens 
the line of sight to the source. Subsequently, at the end of 
the main-on state the inner part of the disk (that may be sur- 
rounded by a hot rarefied corona) covers the source from 



the observer. This picture seems to be confirmed by spectral 
analysis. The opening of the X-ray source is accompanied 
by notable spectral changes including evidence for the pres- 
ence of a strong absorption, whereas the decrease in X-ray in- 
tensity o ccurs more slowly and without appreci able spectral 



changes dG iacconi et al.lll973t iKuster et 
Oosterbroek et alj2000b . 
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2005, see however 



The transition from the off to the main-on state is usually 
called the turn-on of the source. An interesting fact is that 
turn-ons occur near orbital ph ases or b — 0.2 or . 7. Thi s be- 
havior has been explained by iLevine & Jerniganl d 1982b and 



Katz et al.1 dl982l) . These authors assumed that the disk under- 



goes a wobbling with half synodal period due to tidal torques. 
At or b - 0.2 and 0.7 the angle between the disk plane and the 
observer's line of sight changes most rapidly. 

From the data described above it follows that the dura- 
tion of the 35 d cycle is a half-integer number of orbital cy- 
cles. Observations show that in most cases the duration of 35 d 



cycles is 20, 20.5 or 21 orbital cycles dStaubert et al.1 11983 



2006a). 



Another interesting feature of the X-ray light curve is a 
sudden decrease in flux accompanied by significant spectral 
changes - the X-ray dips. The dips were obse rved many times 



by different instruments dShakura et al.lll999l and references 
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therein). X-ray dips can be separated into three groups: pre- 
eclipse dips, which are observed in the first several orbits after 
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Fig. 1. Her X-l ASM light curve sample. 



X-ray turn-on, and march from a position close to the eclipse 
toward earlier orbital phase in successive orbits; anomalous 
dips, which are observed at or b = 0.45 - 0.65 and post-eclipse 
recoveries, which are occasionally observed as a short delay 
(up to a few hours) of the egress from the X-ray eclipse in the 
first orbit after turn-on. 

Several physical models have been proposed to explain the 
origin a nd properties of the X-ray dips. For example, in many 
models dCrosa & Bovntonl ll98ol iGorecki et all 1 19821; llxahv 



19971: ISchandlll 19961) the dips are caused by thickening of the 
outer rim of the accretion disk due to interaction with the accre- 
tion stream; 



are prod uced by the accretion stream ab ove and below the disk 
(see also Bochkarev & Karitskaial[l989l) . However, these mod- 
els made various assumptions about the disk parameters (for 
example, in the original model by Crosa and Boynton the disk 
size was uncomfortably large to explain the X-ray dip march- 
ing behaviour). 

Another model for X-ray dips was proposed by 
1991 In this model pre-eclipse dips occur 



Shakura et al 



when the accretion stream crosses the observer's line of sight 
before entering the disk. This can happen only if the stream 
is non-coplanar to the systems orbital plane. The reason for 
the stream to move out of the orbital plane is non-uniform 
X-ray heating of the optical stars atmosphere by the X-ray 
source, which produces a temperature gradient near the inner 
Lagrange point. The non-uniformity of the heating comes 
from the partial shadowing of the optical star surface by the 
accretion disk. Furthermore, such a stream forces the outer 
parts of the accretion disk to be tilted with respect to the 
orbital plane. The tidal torques cause the disk to precess in the 
direction opposite to the orbital motion. Due to tidal torques 
and the dynamical action of the accretion stream the outer 
parts of the disk develop a notable wobbling (nutational) 
motion twice the synodal period. For an observer, the X-ray 
source can be screened by the outer parts of the disk for some 
time during the first orbit after the X-ray turn-on. This causes 
anomalous dips and post-eclipse recoveries. 



The analysis of Her X-l observations obtained with the 
All Sky Monitor (ASM) on bo ard of the Rossi X-ray Timing 
Explorer (RXTE) satellite (see lBradt et al.lll993l) over the pe- 
riod 1996 February to 2004 September allowed to improve 
the statistics of the averaged X-ray li ght curve of Her X 



1 with respect to previous analyses ( Shakura et alj Il998at 



m 



Scott & Leahy 1999; Still & Boyd 2004). In addition to well- 



known features (pre-eclipse dips and anomalous dips during 
the first orbit after X-ray turn-on) the averaged X-ray light 
curve shows that anomalous dips and post-eclipse recoveries 
are present for two successive orbits after the turn-on in the 
short-on state . These details have already been reported by 



Shakura et alj dl998al) . Now using about two times more data 



we confirm these features to be stable details of the light curve 
(they are also seen in one RXTE/P CA pointed observatio n of 
the source during the short-on state dlnam & Bav kal 2005)). 

In this work we present averaged X-ray light curves of Her 
X-l obtained from the analysis of ASM obser vations. Based 



on this data we modify the model developed bv lShakura et al 



19991) in order to reproduce the details of the averaged light 



Vrtilek & Halpernl d 1985b suggested that the dips 



curve which were left unexplained previously. We show that 
the accretion disk needs to change its tilt during the 35 d cycle 
from ~ 20° in the main-on to ~ 4° in the short-on. 

2. Processing and analysis of the X-ray light curve 

For the analysis of t he X-ray light cur ve of Her X-l we use 
data from the ASM dLevine et a 1. 1996). The archive contains 
X-ray flux measurements in the 2-12 keV band, averaged over 
~90 s. The monitoring began in February 1996 and continues 
up to date. The archive is public and accessible on the Internet 



X-ray data analyzed in this work contain some 50000 
individual flux measurements. They cover ~90 35 d cycles. 
Preliminary processing of t he X-ray light curve wa s carried out 
by the method described bv lShakura et alj dl998al) . 

The goals of this processing are the reduction of the disper- 
sion of the flux through rebinning, resulting in a smoothed light 
curve, and the determination of the turn-on time of each indi- 
vidual 35 day cycle, with the aim to classify the cycles into two 
classes: turn-on around orbital phase ~ 0.2, and turn-on around 
orbital phase ~ 0.7 (they appear with about equal probability) 

The smoothing of the light curves is achieved by combin- 
ing appropriate groups of data points into an aver age (following 
the procedure described in lShakura et alj dl998al) ). The method 
used to determine the turn-on time of the cycles without data 
gaps (we cal l them "good" cycles ) is also similar to the one 
described bv lShakura et al.l dl998al) . For cycles with data gaps 
longer than one day (we call them "bad" cycles) a different 
method is used to determine the turn-on time. Each "good" cy- 
cle was fitted with an analytical function of the following form: 



fit) = A 



1 



1 - exp ■ 



1 +exp((f-A I )/A 2 )\ r A 4 ) 

where t is the time, Ao, .., A\ are the fitting parameters. This 
function consists of two smoothed "steps", the first represent- 
ing the beginning of a main-on state, the second representing 
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Fig. 2. The same sample as in the Fig. [TJ after processing. Vertical 
solid lines mark the time of X-ray turn-ons. 



the end of the main-on. The analytical curves obtained from 
fitting "good" cycles were used as a template library, assuming 
that each "bad" cycle can be well described by one of these ana- 
lytical curves. For a given "bad" cycle we find a corresponding 
curve from our library that fits it best in terms of x 2 , and we 
determine the turn-on time from this curve. 

Figure Q] shows an example of the original X-ray light 
curve. Figure |2] represents the same light curve after process- 
ing. 

2.1. Averaged X-ray Light Curves 

Using the RXTE/ASM (with ~5cts/s from Her X-l in the 
main-on) details cannot be explored in individual 35 day cy- 
cles. However, through superposition of many 35 d light curves 
common details (e.g. X-ray dips, post-eclipse recoveries) be- 
come recognizable. 

This has been done e.g. by Shakura et al.l dl998al) : 



Scott & Leahy! dl999l) : IStill & Bovdl (|2004|). The ASM archive 



has grown considerably, allowing to construct averaged light 
curves with smaller dispersion than in previous works. 

As mentioned above, in most cases turn-ons occur near or- 
bital phases 0.2 and 0.7. Thus all 35 d cycles have been divided 
into two groups - with the turn-on near (p m \, = 0.2 and near 
</>orb = 0.7. Inside each group the light curves were superposed 
and averaged, after shifting them in such a way that the eclipses 
coincided. The superposed light curves are shown in Figs.|3]and 


In the averaged light curve corresponding to the turn-on 
near o ,b = 0.7 (Fig. [4]i one can see that the flux starts to in- 
crease near o ,b = 0.2. But around ol -b = 0.5 this increase 
is suppressed by a strong anomalous dip. This is due to the 
fact that some cycles classified to turn-on near or b = 0.7 have 
some increased flux near or b = 0.2 (just before the anoma- 
lous dip), but statistically insignificant to be recognized in the 
individual light curve. Only when the increase of the flux be- 
fore the anomalous dip is recognizable in the individual light 
curve we consider the turn-on as occurring near or b = 0.2. 



Fig. 5. 3-20 keV PCA light curve of Her X-l during a short-on 
state. Arrows show anomalous dips (AD) which are present on 
two successive orbits after the beginning of the short-on state 

This shows some limitations of determining the turn-on times 
from RXTE/ASM data. 

In the short-on state complicated dip patterns can be ob- 
served. (Fig. |5]l. Anomalous dips near orbital phase or b = 0.5 
are present on two succe ssive orbits after the beginning of the 
short-on state (see also | Leahv et al. 1 120001 lOosterbroek et al 
20001: llnam & Bavkalll2005l) . 



2.2. Statistics of individual 35 d cycles 

The behavior of the duration of the 35 d cycle can be studied 
using the so-called (O - C) diagram - the difference between 
the observed and calculated times (assuming a constant period) 
of the turn-on as a function of time (or 35 d cycle number). 
There is a strong correlation between (O - C) and the pulse 
period development, both on long (years) and sh ort (months) 
time scales (see also Staubert et al. 2000l 2006b ). In addition, 
a lag is observed in the sense that structures in the curves tend 
to appear first in the pulse period curve and then (about 35 days 
later) in the (O-C) diagram. This provides support to the model 
described here (see Discussion). A comprehensive analysis of 
the behavior of the 35 d cycle, pulse period development (P and 
P) and X-ray flux wil l be presented in a forthcoming paper by 
Staub ert et al. I d2006al) . We will show that the correlation to X- 
ray flux is comparatively weak. The statistics of the duration of 
the 35 day cycles is the following: about 40% have a duration 
of ~ 20.5 binary cycles, the rest is equally distributed between 
cycles with ~ 21.0 and ~ 20.0 binary cycles (including a few 
cycles of duration ~ 19.5 binary cycles). 

3. Evidence from the X-ray light curve for a change 
in the tilt of the disk 

The fact that during the short-on state of Her X-l more absorp- 
tion dips appear on the X-ray light curve (the anomalous dips 
and post-eclipse recoveries during two successive orbits of the 
averaged short-on as well as in an individual short-on observed 
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Fig. 3. Averaged X-ray light curves of Her X- 1 corresponding to cycles with turn-ons near orbital phase ~0.2. Vertical lines show intervals of 
X-ray eclipses. 
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Fig. 4. The same as in the Fig.[3]but for cycles with turn-ons near the orbital phase ~0.7. 



byPCA:Fig.|5]and|l nam & Bavkalll2005l) suggests that the an- 
gle e between the disk and the line of sight remains close to 
zero during the first two orbits after the turn-on in the short- 
on state. This is in contrast to the main-on state, where such 
features are observed only during one orbit after turn-on. This 
observation implies a generic asymmetry between the begin- 
ning of the main-on and short-on states. This asymmetry can 
take place only if the disk tilt, 9, changes with precession phase 
(f> pie . Indeed, the angle e of the plane of the disk with the line of 
sight is 

sin e — cos i cos 9 + sin i sin 9 cos pre . ( 1 ) 

where i is the binary inclination angle. So with constant i and 9 
\deld(f> ws ,\ ~ | sin0 pre | and are the same for precession phases 
with e = (i.e. at the beginning of the on-states), although the 
line of sight to the X-ray source lies nearer to the plane of the 
disk during the short-on state. As the binary inclination angle 
; could not change significantly over the precession period, the 
only way to recover the asymmetry is to suppose that 9 is a 
function of d> me . 



The physical reason why this could be so is discussed be- 
low. The disk tilt is generally determined by the action of var- 
ious torques acting on the disk. These include (1) the tidal 
torque from the optical star, (2) the dynamical torque caused 
by gas streams, (3) the viscous torque, (4) the torque caused 
by the radiation pressure and possibly (5) a disk wind torque. 
The joint action of these torques determines the shape of tilted 
and twisted accretion disk (see, for example, calculations of 
the stationary disk shape for the disk wind model in Schandl 
& Meyer 1994). However, if the mass transfer rate through the 
inner Lagrange point into the disk changes with the precession 
phase, so will the tilt of the disk. For example, if free preces- 
sion of the neutro n star is ultima t ely responsible for the 35-day 
cycle in Her X-l ferechedll972tlTru~mper et al.lll986t IShakural 



19951: IShakura et al1ll998bl) . the conditions of X-ray illumina- 



tion of the optical stars atmosphere will periodically change 
with the precession phase. This in turn will lead to changes 
in the velocity components of the gas stream in the vicinity 
of the inner Lagrangian point and hence in the matter supply 
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10 15 
Time (Orbital cycles) 

Fig. 6. The effect of changing tilt of the disk on the angle e 
between the line of sight and the outer plane of a flat disk (solid 
lines) for different phase delays A = 0, 30° and 60°. The phase 
delay causes the asymmetry in the short-on and main-on states. 
The larger A, the more asymmetric is the changing of e. The 
dashed sine-like curve shows the changing of e for the inner 
parts of the twisted disk. 



rate to the accretion disk. Recently, X-ray pulse profiles evo- 
lution with the 35 d phase was successfully reproduced both in 
the main-on and short-on states in the model of freely precess- 
ing neutron star with complex surface m agnetic field structure 
( Ketsaris et al.ll200ot Wilms et al. 2003 ). As a first approxima- 
tion the periodic change of the tilt of the disk with precessional 
phase can be parameterized in the form 



6 = 8o[a + b cos(w pre f + A)] , 



(2) 



where <w pre is the 35 d precessional frequency, A is the phase 
delay, and a and b are numerical constants. The phase delay 
A is zero for the observer whose line of sight coincides with 
the extrema of the disk tilt 9, which is quite improbable. So in 
general one expects A + 0. The effect of the disk's tilt changing 
with the angle e is shown in Fig. [6] for different phase delays 
A = 0, 30° and 60°. It is seen that in a fairly broad range of A 
the possibility appears for e to be close to zero during several 
orbits in the short-on state. 

The anomalous absorption dips and post-eclipse recover- 
ies can be produced by the wobbling motion of the outer disk. 
The wobbling of the outer parts of the disk occurs twice dur- 
ing the orbit al binary period and is primarily caused b y the 
tidal torques ( Katz et al]ll982l iLevine & Jerniganlll982 ). The 
tidal torque vanishes when the companion crosses the node 
line of the disk, so a wobbling appears as a sine-like compo- 
nent with one half of the synodic orbital period on top of the 
mean disk tilt. This explains, in addition to the occurrence of 
the turn-ons at the preferred 0.2 and 0.7 binary phases, the ap- 
pearance of anomalous dips during the first orbit after the turn- 
on both in the main -on and in the short-on states (see Fig. 5 in 
Shakura et alJll999h . 



The appreciable tidal wobbling (with the amplitude of or- 
der of one degree) can be experienced by the disk with the outer 
radius ~ 0.3a (a is the binary separation). Tidal forces also 
cause precession of the disk retrograde to the orbital motion. 
Acting alone, these forces would result in a precessional pe- 
riod somewhat shorter than observed. To explain the observed 
precessional period of the disk (about 20.5P or b, where f or b is 
the orbital period), braking torques must be added. For exam- 
ple, acc ording to our qualita tive accretion stream-disk coupling 
model dShakura et a l] 1 1999b . the motion of the outer parts of 
the disk is driven by the joint action of tidal forces and dy- 
namical pressure induced by the impact of the gas stream on 
the accretion disk. The dynamical action of the gas stream over 
one binary period is opposite to the tidal force so the streams 
can slow down the tidal precession to the observed value. As 
mentioned above, an important feature of our model is that 
the stream is non-complanar with the orbital plane, which is 
due to a complex shape of the shadow produced by the tilted 
and twisted accretion disk on the atmosphere of the optical 
companion (HZ Her). In addition to braking the precessional 
motion of the disk, the stream exerts a torque tilting the disk, 
which counterbalan ces the viscous torque and t herefore limits 
the tilt of the disk (iBardeen & Pettersonlll975l) . So the outer 
regions of the accretion disk are normally kept tilted at some 
equilibrium angle defined by the joint action of viscous and 
dynamical torques (for an alternativ e explanation of the dis k 
tilt induced by the coronal wind see ISchandl & Meverlll994l) . 
Our model has allowed us to successfully reproduce the phase 
behavior of pre-eclipse dips as the result of regular screen- 
ing of the X-ray sou rce by the out-of-the-orbital-plane stream 
JShakuraet al.lll999l) . 



It is important to note that in this model the gas stream im- 
pa cts the tilted disk de eply inside its outer radius (see Fig. 2 
Shakura et al. 19991) at distances r lmp ~ 0.1a from the neu- 



m 



tron star. Assuming the standard accretion disk, the viscous 
time from the impact region to the disk center is around 10 
days, and scales as r 3 ^ 2 . However, the region of the disk be- 
yond n mp , where is no matter supply (called the "the stagna- 
tion zone") and which mediates the angular momentum trans- 
fer outwards for accretion to proceed, can react to perturbations 
induced by the stream impact much faster than on the viscous 
time scale. Indeed, in a binary system, tidal-induced standing 
structure s can appear i n the outer zone of the accretion disk 
(see e.g. lBlondi and references therein) and perturba- 

tions of angular momentum can propagate through this region 
with a velocity close to the sound speed (while the matter will 
accrete on the much slower viscous tim e scale!). Recent analy- 
sis of broad-band variability of SS 433 ( Revnivtsev et al J2005 ) 
suggests that such a picture is realized in the accretion disk in 
that source. 

The physical picture of the change of the disks tilt adopted 
here is as follows. The rate of mass transfer M supplied by HZ 
Her (which for a given moment can be different from the accre- 
tion rate of the neutron star because the viscous time scale for 
mass transport through the disk both delays and smoothes the 
mass flow.) changes periodically over the precession cycle in 
response, for example, to changing illumination of the optical 
star atmosphere due to free precession of the neutron star. The 
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streams action causes the outer disk to precess slower than it 
would do if only tidal torque was acting, and it also changes 
the outer disk's tilt on the dynamical time scale of about 10 
days (see below), which is sufficient to explain periodic disk 
tilt variations over the precession cycle. The wobbling of the 
outer disk is mainly due to tidal forces (the dynamical action 
of the stream provides minor contribution). The viscous time 
scale in the "stagnation zone" is much longer, so the viscous 
torques cannot smooth out the external disk variations. 

In Appendix we present details of the calculation of the disk 
wobbling under the action of tidal and stream torques, and the 
result is shown in Fig. [7] for the standard accretion disk semi- 
thickness H/R = 0.04. In these calculations, there are two char- 
acteristic radii: the outer radius of the disk r out ~ 0.3a which 
determines the tidal wobbling amplitude, and the effective ra- 
dius r e ff ~ 0.18a which determines the mean precession motion 
of the disk. The value of this effective radius is found from the 
requirement that the net precession period of the entire disk be 
equal to the observed value 20.5P or b and the dynamical time 
scale fd = Md/M (which characterizes the dynamical action of 
the stream) be 10 days. This value for the dynamical time is 
chosen to be of order of the viscous time from the impact ra- 
dius ~0. la to meet (quasi) stationary conditions: at larger ?d the 
disk would rapidly empty out, at smaller fd the matter would be 
stored in the disk. We emphasize again that the changing tilt of 
the outer parts of the disk (at r ou t) occurs very rapidly (on the 
sound velocity time scale!), while the precessional motion of 
the entire disk changes on much longer dynamical time scale t& 
in response to the M variations. 

In Fig.|7]we show the result of calculations (see Appendix) 
of the angle e between the direction from the center of the neu- 
tron star to the observer and to the outer parts of the accretion 
disk with non-zero thickness. The complex shape of the disks 
wobbling is clearly seen. The sine-like dashed curve shows 
schematically the angle e for the inner accretion disk regions 
which eclipses the X-ray source at the end of on-states. The 
horizontal line is the observer's plane and the vertical dashed 
lines mark centers of the binary eclipses. The main-on and 
short-on states are indicated. The source is screened by the disk 
when the observer is in the magenta area or between this area 
and the inner disk line. It is seen that the wobbling effects can 
be responsible for the observed (several) anomalous dips and 
post-eclipse recoveries at the beginning of the short-on state. 
These additional dips become more pronounced with the disk's 
tilt decreasing and possibly we see them as unusual anoma- 
lous dips and post-eclipse recoveries at the end of the short-on 
preceeding the anomalo us low state observed by Be ppoSAX 
dOosterbroek et al.ll2000b and RXTE dStill et al.ll20oTb . 



4. Discussion 

For constructing the averaged X-ray light curves of Her X- 1 we 
used the RXTE/ASM data covering ~90 35 d cycles. The large 
amount of data assures good statistics allowing to see some new 
details in the light curves. It is clearly seen that anomalous dips 
and post-eclipse recoveries exist at two successive orbits after 
the turn-on in the short-on state. In our model these features 
are produced by the wobbling of the outer parts of the accre- 



s 

o 



ft 
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Fig. 7. The angle e between the direction from the center of the 
neutron star to the observer and to the outer parts of the accre- 
tion disk with non-zero thickness. The dashed sine-like curve 
shows schematically the angle e for the inner accretion disk re- 
gions which eclipses the X-ray source at the end of on-states. 
The vertical dashed lines mark centers of the binary eclipses. 
Arrows mark anomalous dips (A) and post-eclipse recoveries 
(PE). The binary system inclination i, parameters 6q and A of 
the disk tilt change (Eqf2]i and twisting angle of the disk (de- 
termining the location of the inner disk regions in this figure) 
are chosen such that to reproduce the observed durations of the 
main-on and short-on states and the appearance of two anoma- 
lous dips and post-eclipse recoveries in the beginning of the 
short-on. 



tion disk which cover for some time the central X-ray source 
from the observer. To account for the observed behavior of the 
anomalous dips and post-eclipse recoveries we have to assume 
a change of the tilt of the outer parts of the disk with preces- 
sional phase. In the short-on state the tilt is three times smaller 
than in the main-on state. 

The model we used to calculate the wobbling of the outer 
parts of the accretion disk in Her X-l is highly idealized. 
Although detailed numerical calculations of precessing twisted 
accretion disks includin g viscous forces have been presented 
by se veral authors (e.g. Bardeen & Pettersonl 1975 : Petterson 
1977b . only simplified models have been use d to account for the 
wobb ling motion of the accretion disk (see lLevine & Jernigan 
19821 where only tidal wobbling is considered). Thus, it is 
hardly now possible to check our model, where both tidal 
forces and dynamical actions of the accretion stream are in- 
cluded, by numerical calculations. However, even this sim- 
plified model explains most properties of the averaged X-ray 
light curve of Her X-l: the duration of the main-on and short - 
on states, the appea rance and behavior of pre-eclipse dips 
dShakuraetal.ll 1999b and of anomalous dips and post-eclipse 
recoveries. 

In our calculations we ignored viscous torques exerted on 
the outer regions of the disk. Taking them into account com- 
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plicates the model and can only be done introducing new as- 
sumptions about the shape of the twisted and tilted accretion 
disk in Her X-l. We expect that the effects of the viscosity 
will be the following. The viscosity always tends to reduce the 
tilt of the disk. It counterbalances the dynamical action of the 
stream tending to restore the tilt, so at any given time there is 
an equilibrium tilt angle. Weakening the stream (for example, 
by decreasing the accretion rate M) would ultimately result in a 
smaller tilt of the disk, which we believe is the physical reason 
for the appearance of anomalous low states. 

In the above picture the time needed by the matter to travel 
from the stream impact region r ~ 0.1a to the neutron star (~10 
days) is shorter than the viscous time of the disk at the effective 
radius r e ff ~ 0.18 binary separation (-30 days). So the model 
predicts that the neutron star reacts to the accretion rate varia- 
tions earlier than the precessing accretion disk. This can man- 
ifest itself, for example, in the neutron star period variations 
preceding changes in the accretion disk precessi on period. The 
analys is of long-term correlations in Her X-l (iStaubert et al. 
2006al) indeed suggests such a delay between the turn-on be- 
havior of the source on the O-C diagram and the pulse period 
variations and hence lends support to the model under discus- 
sion. 

5. Conclusions 

The following main results have been obtained from analysing 
and modeling RXTE/ASM X-ray light curves of Her X- 1 . 



1. 



The shape of the averaged X-ray light curves is deter- 
min ed more accurately th an was possible previously (e.g. 



by IShakura et al.lll998al) . This has allowed us to confirm 



the significance and stability of the features in the averaged 
short-on light curve which have been found previously: the 
anomalous dips in two successive orbits after the beginning 
of the short-on state with a post-eclipse recovery after the 
first eclipse. 

2. We argue that the changing of the tilt of the accretion disk 
with 35 d phase is necessary to account for the appearance 
of these features and the observed duration of the main-on 
and short-on states. 

3. We show that the joint action of tidal forces and the dynam- 
ical interaction of the gas stream with the outer parts of the 
accretion disk leading to its wobbling can be responsible 
for the observed behavior of the disk. 
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Appendix A: Wobbling of the outer parts of the 
disk 

The new features discussed in this paper (anomalous dips and 
post-eclipse recoveries in the first two successive orbits in the 
short-on state) call fo r modification of the model proposed by 
Shakura et al. We argue that in order to explain these 

new features, the tilt of the outer regions of the disk should 
regularly change over the 35 d cycle. The physical reason for 
this can be a decrease of the mass transfer rate M supplied by 
the optical component as is expected in our model of a freely 
precessing neutron star as the underlying clock mechanism for 
the 35 d cycle in Her X-l dShakura et al.lll999l) . 



To quantify this behavior, we calculated the net wobbling 
of the disk p roduced by tidal forces a nd the accretion stream. 
In contrast to lSchandl & Meverl (119941) . we ignored the viscous 
torques on the twisted accretion disk (see discussion in Section 
4). 



Calculations were done in the same way as in lShakura et al 



19991) allowing for the change of the tilt of the disk 



To account for the tidal torques, we approximated the 
outer disk as a so lid ring of radius r& (see discussion by 
Shakura et al .119991) . In the quadrupole approximation, the ring 
precesses around the orbital angular momentum with the fre- 
quency 



a> t ~ — ■ 



3 r 3 J 2 cos 6 



(A.l) 



4 V^T? 

where q = M x /M 0.64 is the binary mass ratio and a>b = 
2n/P\, is the orbital binary frequency. 

We denote the precession angle of the disk by <f> (counting 
along the orbital rotation) and the tilt angle of the disk by 9. The 
tidal wobbling of the disk on top of its slow precessional mo- 
tio n in the leading order in u >, / '(2o)h + a> t ) <K 1 can be described 
as ( ILevine & Jerniganlll982h 



d(p 
dt 



1 



2w b 



2wb + <^t 



cos2(0b 



fpre 



d0 
~dt 



— oj t sin ( 



2u h 



sin 2(0b 



^pre 



(A.2) 



(A.3) 



2wb + <^t 

where (f>\, is the orbital phase and pre is the precessional phase. 

Let us now discuss the dynamical action of the stream in 
more detail. The angular momentum of the outer parts of the 
disk Kd are changed by the stream according to the equation 



dK d 

—— = M[r x v] , 

dt 



(A.4) 



where r and v are the distance from center of the disk to the 
impact point of the stream and the gas velocity at the impact 
point, respectively. They depend on the binary phase. In our 
model M is also a function of the orbital phase because the 
matter supply through the inner Lagrangian point depends on 
its position relative to the shadow produced by the tilted and 
twisted disk in the atmosphere of the optical star. 

We will use the right-hand Cartesian reference frame 
rigidly connected to the accretion disk and consider it as non- 
rotating (this assumption is justified by the relatively slow pre- 
cessional motion of the disk). The origin of the frame is at the 




Fig. A.l. The reference frame connected to the disk. The origin 
of the frame is at the center of the neutron star. The vector of 
the angular momentum K of the disk always lies in the z - x 
plane. 

center of the neutron star, the z-axis is normal to the orbital 
plane and the x-axis lies in the orbital plane and is directed 
such that the disks angular momentum vector always lies in the 
Z — x plane (see Fig. IA. Il l: 



K d = (/ST d sin6»,0,/i:dCos6»). 



(A.5) 



Then the change of the disk momentum c/Kj = 
(c/K x , c/K y , dK z ) can be written through the change in its 
components as 



dK d = dK x sin + dK z cos 6 , 
K A 6 = dK z cos 8 - dK z sin 9 , 
K,\ sin 9d<p = dK y . 



(A.6) 



The time evolution of the disk angular momentum components 
(Eq. (lA.4t ) in this frame now can be written in the form: 



dK x 
dt 
dK y 

~~dT 



Mv K \a(yv z - zv y ) , 



Mv K ia(zv x - xv z ) , 



(A.7) 



dK z 

= Mv le ia(xvy - yv x ) , 



where x, y, z are dimensionless coordinates of the stream im- 
pact point (in units of binary separation a) and v x ,v y ,v z are 
the components of the velocity of the stream at the impact 
point (in units of relative velocity of binary components v re i). 
For computing ballistic trajectories of the gas stream from 
the Lagrangian point we used the primed right-hand reference 
frame x',y',z' rigidly rotating with the binary system. The ori- 
gin of the frame is at the center of the neutron star, the z'-axis 
is normal to the orbital plane and the x'-axis is pointed to the 
center of the normal component. The transition from that frame 
to the one of the precessing disk is made by usual coordinate 
transformation: 
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x — x' cos <5 - y' sin <E> , 
y = x sin d> + y' cos Q> , (A.8) 
Z = z', 

and for velocities: 

v x = v' x cos O - v' y sin <t> - J(x') 2 + (y') 2 sin <1) , 



v y = v x sin O + v y cos <1> + ^J(x') 2 + (y') 2 cos <1> , (A.9) 

where <t> is the angle between the coordinate axes x and x' . 

It is convenient to rewrite all equations in a dimensionless 
form. The angular momentum of the disk can be written as 

|K d | = yM d w k 4, (A. 10) 

where r d is the effective disk radius (see Section 3), M d is 
the mass of the disk and co k = jGM x /r^ s is the Keplerian 
frequency at the effective radius and y is a numerical co- 
efficient accounting for the surface density distribution 
for example, in the standa rd Shakura-Sunyaev accretion disk 
dShakura & Sunyaevll 19731) S oc r 3/4 and y = 5/7. 



Now, from equations ( IA.6t and ( 1A.71 i) and making use of 
the Kepler's 3rd law, we arrive at the following equations for 
the change of angles (f> and 9 with the synodic angle r = tu s 

(iO s = COpre + W b ): 



d<f> M \M X +M Q fa 

sm6— = — — J — — J — [zv x -xv z ] , 

dr yM d w s \ M x \ r d 



</0 M M x + M a nr 

T = a/ 77 A / — [cos 6»(yv z -zv y )- (A.ll) 

- sin ^(xvy - yv x )] . 

Using these equations, we calculated the change in angles <p 
and (A<f> and A6) step by step and sum them up to find their 
behavior over one precessional cycle. The initial components of 
the velocities v' x , v' , v' z were chosen as a f unction of the s ynodi c 



phase O in the way described in detail bv lShakura et al.l (11999b . 

The free parameters are the mass of the accretion disk M d , 
the maximal disk tilt do and its change, the disk semi-thickness 
/z/r d at the outer boundary and the binary inclination angle 
;. We fixed the angle i = 88?6 and the disks semi-thickness 
h/r<i = 0.04 (corresponding to a disk opening angle of about 
5°). The disks tilt was allowed to change with precessional 
phase according to the Eq. (fJJ. 

The angle 9q = 20° and numerical coefficients (a = 0.6, 
b = 0.4 and A = 15°) in Eq. (f2]i were chosen to reproduce the 
observations. Now we can calculate the angle e between the 
line of sight and the outer parts of the disk making use of Eq. 
0- 

In contrast to the case with constant tilt considered ear- 



lier by IShakura et al.l (119991) . the change of e with precessional 
phase is determined also by the change of the tilt 6((p pie ) [Eq. 
10]. The effect of the phase delay A is illustrated by Fig. [6] 
the larger A is, the more asymmetric is the changing of e The 
results are shown in Fig. [7] 



